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SUMMARY 


A  number  of  estimates  and  tests  for  mean  life  and  other 
parameters  derived  under  the  exponential  distribution  assumption 
are  studied  under  the  alternative  condition  that  the  distributio  .  has 
increasing  (decreasing)  failure  rate.  The  estimates  considered 
are,  for  the  most  part,  based  on  censored  and  truncated  samples. 
It  is  shown  that  these  estimates  generally  favor  the  producer  (con¬ 
sumer)  in  the  IFR  (DFR)  case.  Properties  of  order  statistics  and 
their  spacings  from  distributions  with  increasing  (decreasing) 
failure  rate  are  presented. 
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1.  Introduction.  In  f  ii  nda  nu- 1 .  t  u  1  paper  in  »h«-  literature  of  life 
testing  Fpstein  ano  Sohel  (IdSI)  introduce  li^^e  tc-st  procedu/es  based 
on  the  exponcT.t  iri  1  distribution.  Ihc-se  procedures  have  been  codified 
in  a  Department  of  nefense  nandbook  (1960)  and  arc*  nou  idely  em¬ 
ployed.  7elen  and  Dannemiller  (l'<61)  shoo  bv  sampling  from  Weibull 
distribution  alternatives  that  these  procedures  are  not  robust  in 
testine  for  mean  life.  Hovvever,  as  Antelman  and  Savage  (unpub.  ) 
have  pointed  out  they  mav  be  robust  in  testing  for  cert..n  percentiles. 
For  certain  loss  functions  based  on  percentiles,  these  procedures 
seem  to  be  robust.  Since  statistical  procedures  bas'*d  on  the  ex¬ 
ponential  distribution  have  a  great  deal  of  intuitive  appeal  and  com¬ 
putational  simplu  itv  we  investigate  their  properties  relative  to 
alternative  distributions  having  incr«*asing  failure  rate  (IFR)  or  de¬ 
creasing  failure  rate  (DFR). 

This  paper  essentiallv  confirms,  theoretically  and  more 
general!'',  the  sampling  results  of  the  7  c*  len  -  Da  nne  mi  lie  r  paper  for 
statistics  derived  under  the  exponential  assumption.  Using  Weibull 
distribution  alternatives  (vcith  parameter  values  vc-hich  insure  that 
the  distribution  has  increasing  failure  rate)  7  elen- Dannemiller  show 
that  the  use  of  these  statistics  mav  result  in  substantially  increasing 
the  probability  of  accepting  items  having  poor  mean  lives.  We  shov 
that  these  estimates  for  the  mean  are  positively  (negatively)  biased 
when  the  distribution  is  IFR  (DFR).  Also  \>  e  obtain  bounds  on  the 
expected  values  of  the  e>Donential  estimates  ‘or  the  distribution 
function  and  bounds  on  the*  expected  values  of  the  order  statistics. 


In  the  last  sec  tion  varicjus  properties  of  IFR  (DFK)  carder  statis¬ 
tics  are  presentee*. 


Pr elinaina r ie s .  Let  X  denote  a  random  variable  with  righ*  continu¬ 
ous  distribution  F  such  that  F(0  )  =  0.  If  F  has  density  f  then 


r(t)  = 


f(t) 

rT-.Tft)r 


is  known  as  the  failure  rate.  Note  that  r(t)  = 


-  logf  1  -  F(t)l  when  a  density  exists.  For  this  reason,  v  e  say 
that  F  is  IFR  (DFR)  for  increasing  (decreasing)  failure  rate  if 
logf  I  -  F(t)l  is  concave  where  finite  (convex  on  f  0,  or)  ).  Note 

that  any  IFR  (DFR)  distribution  with  specified  mean  can  be  ex¬ 
pressed  as  the  limit  of  continuous  IFR  (DFR)  distributions  '<  ith 
the  same  mean.  Hence  for  many  of  our  results  it  is  sufficient  to  con^ 
fine  attention  to  continuous  IFR  (DFR)  distribution?. 

We  often  use  the  well  known  fact  that  if  F  is  continuous, 
then  Y  =  -  0  logp  (X)  is  exponentially  distributed  with  mean  h 
where  F'(x)  -  1  -  F(x)  .  Repeatedly  we  use  the  fact  that  if  F  is 
IFR  with  mean  0  then  there  exists  Xq  >  0  such  that 


y  =  -0  log  r  (x) 


^  <  X  for  X  <  Xq 
^  >  X  for  x2  Xq 


This  is  evident  from  log  concavity  and  the  bounds  on  IFR  distri¬ 
butions  given  in  Barlow  and  Marshall  (1964).  The  inequalities  are 
rever8»*d  when  F  is  DFR. 

Unless  otherwise  indicated  we  denote  ordered  observations 


from  a  random  sample  of  size  n  based  on  a  random  variable  X 


by  X,  X  .  We  define  X»  -  0  . 

1  —  —  n  0 


1.  Estimates  based  on  censored  samples.  Assume  n 

put  on  life  test  and  let  X,  <  X  ,  <  .  <  X  denote 

i  —  c  —  —  n 

observations.  If  F  has  density  f  such  that 

X 


f(x;  6) 


X  >  0 


‘  0 


X  <  0 


Items  a  re 
the  ord'Tcd 


then 

r  r 

y  X.  Mn  -  r)X^  (n  -  i  f  1)(X.  -  X._,) 

(2.1)  -  "r 

r,  n  ‘  r  r 

(1  <  r  <  n)  is  the  maximum  likelihood  and  minimum  variance  un¬ 
biased  estimate  for  B  based  on  the  first  r  order  statistics 
(Epstein  and  Sobel,  19‘^3). 

The  normalized  spacings  D.  =  (n  -  i  +  1)(X.  -  X.  j)  which  enter 

A 

into  the  computation  of  ^  have  a  natural  intuitive  appeal.  Thev 

have  also  been  used  as  the  basif  for  a  statistic  to  test  for  IFR 
(Proschan  and  Pyke,  in  preparation).  We  shall  derive  and  use  sev 
eral  p.-operties  of  these  spacings  when  F  is  IFR  (DFR)  .  Sinc»- 
the  normalized  spacings  are  independent  and  identically  distributed 
in  the  exponential  case  (Epstein  and  Sobel,  1953),  Theorem  2.  1  belin> 
is  quite  intuitive. 

A  random  variable  X  is  said  to  be  stochastically  smaller 
than  a  random  variable  Y  if  and  only  if  Pf  X  >  x]  <  P[  Y  >  xl  for 


all  X. 
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Theorem  Z .  I.  If  F  is  IFR  (DFR)  the  normalized  spacings 

(n  -  i  +l)(X^  -  j)  are  stochastically  decreasing  (increasing)  in  i  . 

Proof.  Assume  F  is  IFR  and  let  F  (x)  =  1  -  F(x).  Note  that 

PfnX,>  xl  =  (r(-^)r  > 

I 

since  [F(t)]^  is  decreasing  in  t  Let 

F  (x)  =  x)  -  F(u) 


and  note  F^(x)  >  F(x)  Given  that  X^  =  u  is  observed,  X^  -  Xj 
is  distributed  as  the  first  order  statistic  from  a  sample  of  size 
n  -  1  each  with  distribution  F^(x)  .  Hence 


Pf  (n  -  DIX^  -  Xj)  >  X  I  Xj  =  u]  =  f  . 

Conditioning  on  X^  we  have 


P(nX,>x|  =  (F(i)l''> 

=  Pf  (n  -  1)(X^  -  Xj)  >  X  I  X  j  =  ul 
for  all  u  >  0.  Unconditioning 

P[nX^>x]  >  [  [  dG(u) 


=  P[(n-l)(X,  .Xj)>  xl 

where  G(u)  =  1  -  fF(u)l^  is  the  distribution  of  Xj  Hence  we 
have  show'n  that  nX^  is  stochastically  larger  than  (n  -  -  Xj)  . 

In  a  similar  manner  we  can  show  that  (n  -  i  +  1)(X.  -  is 

stochastically  larger  than  (n  -  i)(X^^j  -  X.)  for  i  =  Z,  3 . n 

All  inequalities  are  reversed  for  DFR  distributions.  || 

As  an  immediate  consequence  of  Theorem  Z.  1  we  have  that 
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E/f(n  -  i  +  1)(X,  -  X  I 

I  4  l-l  j 

is  decreasing  (increasing)  in  i  for  ^  >  0  when  F  is  IFR 

A 

(DFR)  Using  this  fact  we  can  show  that  6^  ^  is  positively 
biased  when  F  is  IFR. 

Corollary  2 .  <?.  If  F  is  IFR  with  mean  6  ,  then 

0  <  E[  6  ]  <  n--  lor  r  =  1,  2,  ....  n. 

—  r,  n  —  r 

All  inequalities  are  sharp. 

Proof.  From  Barlow  and  Proschan  (1964a,  p.33)  we  know  that 

E[  6,  ]  >  6  .  Also 

1,  n  —  r 

h(r)  =^|Ef  (n  -  i  +  1)(X.  -  X,  j)l  -  oj 
i=l 

exhibits  at  most  one  sign  change  as  a  function  of  r  since 
E((n  .  i  +  1)(X.  -  X._j)] 

is  decreasing  in  i  by  Theorem  2.1.  But  h(l)  >  0  and  h(n)  =  0  , 

which  implies  h(r)  >  0  for  r  =  1,  2,  .  .  .  ,  n  .  Hence 

E[©  1  >  f)  . 

r,  n’  — 

Clearly  the  bound  is  attained  by  the  exponential  distribution  so  that 
it  is  sharp. 

To  show  the  upper  bound  we  note 
r  n 

V  V 

^  X.  +  (n  -  r)  X  <  ;  X. 

1  r  —  ,  i 

i=l  i=l 

for  every  sample  realization.  Hence 

Ef  r  6  ]  <  n  0 

r,  n'  - 

or  E[0  1  <  n—  . 

'■  r ,  n  '  —  r 


6 


Since  equality  ie  attained  with  distributions  degenerate  at  0  (vhich 
is  the  limit  of  IFR  distributions)  the  bound  is  sharp.!] 

Corollary  Z.  3.  If  F  is  DFR  with  mean  0  ,  then 

0  <  El  G  ]  <  G  for  1  <  r  <  n  . 

-  '  r .  n '  -  - 

All  inequalities  are  sharp. 


Proof.  The  upper  bound  follows  from  Theorem  Z.  1  and  the  method 
of  proof  in  Corollary  Z.Z.  To  show  that  the  lower  bound  is  sharp,  let 


r(x)  =  <j 


r  0 


cx 


€  e 


X  <  0 


X  >  0 


where  e  >  0  is  arbitrary.  Then  F  is  DFR  with  mean  0  and 

i-1 

P[  X.  >  x]  f  Fix)]-’  f  rix)]''"-’ 

j=0 

i-1 


K\  ^  *  j  ) 


J=0 


Hence 


Et  X.) 


fpfx  >x!  dx  .  y 


J  /  (P'j) 


when  0  <  e<  1  Since  e  is  arbitrary  we  see  that 

Ef  X  J  >  0  (I  <  i  <  n) 

is  sharp.  j| 

.th 


For  convenience  we  now  denote  the  i  order  statistic  from 


a  sample  of  size  n  by  X. 

i,  n 


Theorem  Z.  4.  If  F  is  IFR  (DFR)  (n  -  i  +  l)(X.  -  X.  ,  )  is 

-  I .  n  1  -  i ,  n 


stochastically  increasing  (decreasing)  in  n  for  fixed  i.  Hence 

]  .  (1  <  r  <  n) 


(» 


7 


Proof.  Assume  F  is  IFR  .  Let  G.  (x)  =  P[\.  <  xl  and 

-  i,  n  ^  1 .  n  -  ' 


note  that  G.  (x)  <  G.  ,(x)  for  samples  from  any  distributi 
i,  n  -  1,  nfl  ^  ^ 

Now 


1  on . 


Pf  (n  -  i)(X.  ,  -  X.  )  > 

‘  '  1  +  1,  n  i,  n  - 


^  n<“' 

jQ  uri”!  i,n 


-  'q  '  u  n  +  l-i  i.  n 


(u) 


r 


cr 


dC.  ^,(.1) 
1 ,  n  +  1 


1^ 

the  first  inequality  holds  since  f  F(t)l  *  is  decreasing  in  t  when 
F  is  IFR  . 

All  inequalities  are  reversed  when  F  is  DFR.  || 

A 

Thus  when  F  is  IFR  ,  the  estimate,  ,  of  mean  life 

r ,  n 

based  on  a  sample  censored  on  the  right  becomes  worse  with  in¬ 
creasing  n  when  r  (1  <  r  <  n)  remains  fixed. 


Acceptance  Sampling.  Statistical  methods  for  testing  hypotheses 
about  the  mean  of  an  exponential  distribution  depend  on  »he  statistic, 
^  .in  the  case  of  censored  samples  (Epstein,  1960a).  For  test¬ 
ing  the  hypothesis  Hq:  ^  against  the  alternative 

6  =  subject  to  P[  reject  ^  j  true]  =  n  .  the  rejec  - 

tion  region  is  of  the  form 


5  ^ 

r ,  n  - 


2  r 


If  X,  (2r)  <  2r  ,  then  we  shall  prove,  using  Lemma  2.  b,  that 


1 


H 


(2.2) 


P  P  <  -  -- i  ^ -  F  IFR,  0  0^ 

r ,  n  —  Z  r  0 


—  r ,  n  — 


r 

\- 

TT 


F  ( X )  -  1  -  e  ^ 


0 


r> 


so  that  the  exponential  test  provides  a  size  q  test  whf-n  the  failure 
distribution  is  IFR  .  To  see  (Z.Z)  we  need  the  fnllo"inp  easi!v 
verifiet'  result  v'hirh  '•  e  present  without  proof 


Lemma  Z.S.  If  (f)  is  concave  '^(0)  r  0 


X-  ^  0  f i  =  1 ,  Z . n)  ,  then 


and  a.  ■>  1, 
1  — 


n 

^  V 

<ti  ax 

1  1  — 

i=l  '  i=l 


n 

<  ^  a .  <t>(x. ) 

I  I 


Let  <i)  \y)  =-f)log  F(y)  so  that  X.  =  4)(Y.) 

where  Y.  is  the  i^^  order  statistic  of  an  exponentially  distribu 

ted  random  variable  with  mean  f>  .  Then 

r 

A  ^ 

rO  =  >  <t)(Y.)  +  (n  -  r)  6  (Y  ) 

r ,  n  !  r 

irl 


>  <t> 


)  Y.  ^  (n  -  r)Y  ! 


/  M 

1=1 


by  Lemma  Z.S.  Using  the  bounds  on  IFR  distributions  (Barb 


and  Marsnall,  1964)  and  letting  Z  =  F  Y.  t-  (n  -  r)Y  .  v'e  have 

*’  i=l  '  ’’ 

Pp  f<t>(7.*)  <  c]  <  P^f  7*  -  c] 


We  obtain  (2.2)  by  letting  ^  ~ 


when  c  <  0 


>) 


c 


v;  (^r) 
1  -  O' 

— 7 - 


Samplit.p  with  I^rplac  rnicnt .  Suppose  no'*.  that  failed  it'-nis  are  re¬ 
placed  at  failure.  In  this  case  the  bias  of  the  usual  estimate  for  (' 

is  even  greater  than  in  the  non  -  replac  ement  case. 

,  .  t  h  . 

Let  X.  denote  the  time  of  the  i  failure  -Alien  failed  items 


are  replaced.  The  maximum  likelihood  estimate  for  based  on 

tr.e  exponential  assumption  is,  in  this  case 


A  1 

-  f  nX, 
r ,  n  r  1 


•  n(X, 


X,)  .  , 


‘  "'X  -  X-i" 


nX 


r 


(Epstein,  1960b). 

Theorem  Z,  If  F  is  IFR  with  mean  6  ,  then 

f<  ^  K\f>  I  -  FfO*  1  -  Ff  X  1  1  ■'  r  -  n 

—  r,n—  r.n—  r  r  —  — 


Proof.  Clearly  X^  X^  for  any  distribution  F  so  that  the 
upper  bound  is  obvious.  To  sho"-  the  remaininp  inequality  we  intro¬ 
duce  the  following  fictitious  replacement  policy: 

Policy  A:  Replace  a  failed  item  with  a  good  item  of  the  same  "age" 
**  .th 

Let  X.  denote  the  time  of  the  i  failure  under  this  policy  and 

r 

A,  A  *  1  A  A 

1  )  (n  -  1  .  1)(X.  -  X.  ,  )  . 

r , n  r  i  i  - 1 

1=1 


It  is  clear  that  F  f 1  -c  F  ( 1 

r ,  n  —  r.n 

since  under  the  IFR  assumption  the  conditional  mean  life  of  an  aged 
item  18  less  than  the  mean  life  of  a  new  item. 
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We  need  only  sho'*'  F  1  Ff  ^  1  .  Let 

r ,  n  -  r  ,  n 

F  (x|  =  .ind  C(u)  =  P[X,  -  111  .  Then 

“  F(u)  ‘  - 

Ef(n-l)(X,-X)1  T  ^  fF'(  — ^  )  1  ^  dx  dG(u )  Simila  r  Iv 

^  I  ^0  -0  " 


under  Policy  .  . 


pfx""  -  x;"  >  ^  1  X,  =ui  -  fr  (^)i" 

^  1  —  n  '  1  u  n 


and 


Ffn{X^^  -  Xj  )] 


0  >^0 


f  r  (  -  )!"  dx  dG(u)  . 
u  n 


Sine  e 

have 

To  sho\«.' 


f  r  (  — )  1  ^  f  F*  (  — ^  )  1  ^  ^  \*  hen  F  is  IFR  v.e 

u  n  —  ‘  u  n  - 1  ' 


Ff(n  -  1)(X^  -  X^)]  <  Efn(X^’'  -  Xp)l 


r  ,  r  Jjt  )}t  3i|t 

F  (n  -  i  f  l)(X.  -  X.  ,)1  F  (n  -  i  +  1)(X.  -  X. 

i  i- 1  —  i  1-1 


for  2  <  i  <  n  we  proceed  as  above  except  that  the  definition  of  G 
is  different  for  the  two  policies.  For  example,  for  i  =  3  let 

G{u^  I  Uj)  =  Pf  X^  <  I  Xj  =  Ujl 


n  - 1 


and 


=  1  -  f  F  -  u,) 

G  (u^  I  u^)  =  Pf  X^  <  !  Xj  =  Ujl 


=  ‘  ■  f  ■  “l”"  ■ 


Since  1  ''!>  i  G  (u^  I  Uj)  we  have 


II 


Ff  (n  -  ^)(X,  -  X  ,)| 

}  *_ 

•  or  . »  or 

=  \  \  [  F'  (  ■%  )  1  "^flx  dG(M  J  u  )  dG(u  ) 

.  0  .  0  ‘ 

»  or  .or 

^  \  \  (— "^dx  dG  (ii  J  u  )  dG(u.) 

•  or  I  or 

^  \  \  [F  (-iTdx  dG’'(u,  i  u  )  dG(u  ) 

•0-0  ^  ^ 

r  *  ^  ^ 

-  F(  n(X^  -  X,)l  .  A  similAr  ccirdition  ng  arj’Limcnt 
wnrksfor  3''i^n  .|| 

Bounds  on  estimates  for  the  reliability  function.  The  minimum 
variance  unbiased  estiinate  fo  R(t)  =  Ftt)  (t  is  fixed)  under  the 
exponential  assumption  is 

R^(t)  =  maxf  0,(1--^)*^ 

r 

r 

where  7  =  ^  X.  +  (n-r)X  For  a  discussion  of  such  minimum 

1=1  ' 

variance  unbiased  estimates  see  Tate  (19‘i9).  For  convenience, 
assume  that  =  1  Then,  under  the  exponential  assumption,  7  ^ 
has  density 


Theorem  2.6.  If  F  is  IFR  with  mean  6  =  1  and  t  ^  R  =  1, 
then 

Ef  Rj(t)l  >  \  f  1  -  K^fy)  dy  +  (l-t‘/’  f?^(y)dy  . 

Proof.  W  thout  loss  of  generality  we  may  assume  F  continuous 
(see  preliminaries).  Let  (t>  \y)  =  -log  F(y)  Then  ct)  is  con¬ 
cave,  increasing  and  <i>{0)  =  0  .  If  Y.  is  the  i 


order  statistic 


from  an  exponentially  distributed  random  variable  with  0  -  \ 
then  X.  -  <h(Y.)  is  the  i*^  order  statistic  from  an  IFR  random 

1  i 

variable  with  distribution  F  and  mean  P  -  1.  Furthermore 


7 


(1>(Y.)  Mn  -  r)  (f(Y  ) 
i  r 


i  1 


Y.  f  (n  -  r)Y 
1  r 


,  A 

.  f  7.^1 


by  the  previous  lemma.  Therefore 

Rj(t)  -  maxf  0.(1  -  )*" 


ma X  f  0,(1 


■V)'''  1 


<^(Z^) 


Since 


<t>  (y) 


V  -  1 
y  1 


•ve  have 


F-f  R,(t)]  ''  \  f  1  -  7,  1  g_(y)  dy  +  (  (1 

1  -  t  ^  •  1 


gj.(y)  dy 


The  maximum  likelihood  estimate  for  R(t)  under  the  exponen¬ 


tial  assuiiaption  is 


t 


^  ,  s  '  r .  n 

R^(t)  =  e 


where  B  was  defined  in  (2.  1)  .  Pugh  (1063)  has  shov  n  that 

r,  n 

A 

under  the  exponentia.  assumption  **  negatively  biased  when 

tne  true  reliability  R(t)  >  1  -  .  368  .  Assuming  F  is  IFR 


r  ^  ^ 

we  can  obtain  a  lower  bound  on  E|R^(t}! 


Theorem  2.7.  If  F  is  IFR 


with  mean  6  -  1  ,  then 


I  f 


tr 

fv^ 


t  r 

TT 


0 


♦*  Rj.(y)  rfy  ^  \ 


Kj.(v)  dv 


The  proof  parnllels  ihtt  of  Th(‘f)retTi  ^  .  f) . 


1  ^ 


yistimates  and  confidfnu  e  hounds  on  p«  r(  entiles.  If  F 

th 


with  mean  f>  and  p  perc  entile 


then 


is  \VR 


[  -logll  -  p)lf‘  1  ^  p  1  0  ■ 

(See  Barlow  and  Marshall  (196-1)  ).  Hence  bv  Theorem  d.  1 


F  .  0 


f  -l')g(l  -  p)] 


r ,  n 


.  >  f 

I  -  'P 


for 


=  1 ,  Z  ,  .  .  .  ,  n 


■hile 


F  -  6^  ^  f  -logd  -  p)]  '  5  "'p 


and  one  might  be  tempted  to  use  these  estimates  to  bracket 

Intuitively,  we  want  a  confidence  interval  to  ha'  e  small  expect¬ 
ed  width  when  it  covers  the  true  percentile.  The  usual  distribution- 
free  confidence  intervals  based  on  order  statistics  have  smaller 
conditional  expected  width  under  the  IFR  (DFR)  assumption  than 
under  the  exponential  assumption,  given  that  the  interval  contains 
the  true  percentile.  To  see  this  let  Y  =  -  ^p)  ^rid  note 

that  Y  is  exponentially  distributed  with  p*^  percentile  when 


F  is  continuous.  Suppose  that  X.  <  C  <  X 

i  -  p  -  1 


Then  clearly 


Ef  X.  -  X.  j  X.  <  q  ''  X.l  <  Ef  Y  -  Y.  Y.  <  ?  <  Y  ]  . 

j  ili-p-j’-  j  i  i-p_j' 


1> 


3.  Estimates  based  on  truncated  samples.  If  n  items  are  placed 
on  life  test  and  if  sampling  is  terminated  at  time  T  ,  the  associa¬ 
ted  sample  is  called  a  truncated  sample.  I,  et  X ,  X  ,  X 

I  -  c  —  —  n 

denote  an  ordered  sample  from  a  distribution  P'  and  let 

r 

V(T  )  ^  X^  f  (n  -  r)T 
1-1 

where  r  is  a  random  variable  and  denotes  the  number  of  X's  less 
than  T  .  Then  V{T)  isthe  total  life  observed  up  to  time  d' 
This  statistic  occurs,  for  example,  in  sequential  life  tests  for  the 
exponential  case  (Epstein  and  Sobel,  19^^^).  It  is  not  surprising  that 
this  statistic  also  has  greater  expected  value  under  the  IF'R  as¬ 
sumption  . 

X 

■  "B 

For  convenience,  let  G(x)  =  1  -  e  ^  . 

Theorem  3.  I.  If  F  is  IFR  (DFR)  with  mean  0  .  then 
Ep,fV(T)l  >  E^fV(T)l 

Proof.  Assume  F  IFR  and  let  X,  <  X,  <  ....  X  denote  an 
ordered  sample  from  F  Without  loss  of  generality  we  may  as¬ 
sume  F  continuous.  Let  y  =  -fdog  F(x)  We  kno\>  there  exists 
Xq  >  ft  such  that  x  >  -ftlogF(x)  for  x  <  x^  and  x  <  -ftlogF(x) 
for  X  >  Xq  (Barlow'  and  Marshall,  1964). 

Let  Y.  =  -ftlogF(X.)  If  T  X-  .  then 

1  i  —  0 

XV  s 

^X.  f(n-  r)T  >  )  Y.  (n  -  r)T  Y.  f  (n  -  s)d 

—  f  \  —  _  \ 

i  =  l  ml  1  =  1 

where  r  (s)  denotes  the  number  of  X's  (Y’s)  less  than  T. 


1^ 


Hence  for 


i  ’‘o 


I  ,et 


For  T  >  X , 


Fp.[  V(T)1  >  F^f  V(T)] 


Y. 


r 

V 

i=l 

r 


Y 

I  ‘  i 


if  Y  T 
1  — 


T  othe  r  •’  i  se 


s 


i=l 


X.  Mn  -  r)T  -  '  Y.  -  (n  -  s)T 

i  i 

i  =  l 


r 

V  * 

X.  4-  (n  -  r)T  -  '  Y.  -  (n  -  r)T 


i  =  l 


n 

V 


>  )  (X.  -  Y.)  +  y  (X.  -  Y.) 

—  1  1  __  I  I 

i  =1  i =r  f 1 


since  X.  <  Y.  for  i  >  r  ,  Hence 

i  —  1 


F„f  V(T)]  -  E_f  V(T)1  >  E|  /  X.  '  -  f:  ^  Y.  =  0 

r  U  —  1  i 


i=l 


i=l 


for  T  >  X 


0  • 


A  similar  argument  holds  for  the  DFR  case. 
Consider  the  estimate 


;^(T)  =  xiii  = 


X.  +  (n  -  r)T 
1=1  ' 


nT 


ii  r  >  0 

if  r  =  0 


When  F  is  the  exponential  distribution,  ^(T)  is  the  max'mum 
likelihood  estimate  of  .  In  this  case 


E[^(T)l  =  f>  -  ’  > 

1  - exp(- 
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since  r  and  ^(T)  are  negatively  correlated  (Bartholomew, 
1957).  In  the  IFR  case,  this  statistic  exhibits  even  greater  bias 

X 

for  T  0  .  As  before  let  G(x)  =  1  -  e  71  . 

Theorem  3.  2.  If  F  is  IFR  •••ith  mean  6  ,  but  not  degenerate 
then 

Fp,f  6(T)  I  r  >  1]  >  F^[R(T)  |  r  >  l]  for  T  ^  6  . 

Proof.  Assume  F  is  continuous  and  let  Y.  "  -61ng  F(X.)  as 


before  and  let 


b.  = 

i 


r 


if  X.  <  T 

1  — 

otherwise 


if  Y  <  r 
1  — 


We  can  write 


0  other\«isc  . 

n 


T  -  a  (T  -  X.) 


fi(T) 


i=l 


n 

^  - 

/  1 

i  =  l 


if 


n 

\ 

iTi 


a .  ">1 

i  — 


Assume  T  <  P  .  As  in  the  previous  proof  X.  <  T  implies 

Y.  <  T  and  hence  a  <  b.  .  If  a.  =  1  ,  then  b.  =  1  and 

i  i  —  i  i  i 


T  -  X.  <  T  -  Y.  .  Hence  if  F  a.  >  1, 


‘  irl 

n 

Y  r 

)  '  T  -  a.(T  -  X.) 


1  — 


\  ' 


1 


T  -  b,(T  -  Y.) 


i=l 


>  i=l 


n 

L 

1=1 


Ep,  ff>(T) 


n 

a. 

)  b. 

i 

/  i 

1 

i=l 

A 

r  >  1]  >  E„ff5(T) 

r  >  1]  for  T  <  0  . 

—  —  o 

and 
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Inverse  Binomial  Sampling. 

Nadler  {I960)  has  considered  the  following  type  of  sampling. 

An  item  having  life  distribution  F  with  mean  0  is  put  on  test  until 
it  fails  or  time  t  has  elapsed;  at  this  time  the  item  is  replaced  by  a 
fresh  item.  This  is  repeated  sequentially  until  r  actual  failures 
are  observed.  The  number  N  .  of  items  that  have  to  be  tested  until 
the  r  actual  failures  are  obtained  is  a  random  variable.  Nadler 

X 

(1960)  showed  that  when  F(x)  =  1  -  e  ^  ,  an  unbiased  estimate  of 
0  is 

r 

A  1  ^  \  ^ 

(t)  =  -  -  ^  Y.  f  (N-r)t  - 
r'  r  ^  1  '  ' 

i=l 


where  the  Y^,  .  .  .  ,  Y^  are  the  r  life  lengths  not  exceeding  t  . 

We  show  next  that  when  F  is  IFR  (DFR)  with  mean  0  ,  then 

A 

6^(t)  is  biased  high  (low). 

Theorem  3.3.  If  F  is  IFR  (DFR)  with  mean  9  ,  ,then  E6^(t)  >(<)0  . 

Proof.  Let  F  be  IFR  .  Let  Z.  denote  test  time  elapsed  between 

the  i  -  1^^  failure  time  and  the  i^^  failure  time,  i  =  4>  •  •  •  >  r, 

where  the  0^^  failure  time  is  defined  to  be  0  .  Then 

r 


0^(t) 


.  1  V 
■  r  ^ 

i^i 


z.  . 

I 


Next  consider  an  alternate  testing  procedure  differing  in  that 

replacement  occurs  only  upon  failure.  Let  Z/.  =  test  time  elapsed 

1  etween  the  i  -  1^^  failure  and  the  i^^  failure  under  the  alternate 

testing  procedure.  Now  since  F  is  IFR  ,  is  stochastically 

larger  than  Z.  .  It  follows  that 
*  1 


1^- 


^  1  \  I  \ 

Kf)  (t)  =  -  ^  F.7>-^K7.^e 

r  r  ,  1  -  r  i 


The  inequality  is  reversed  -vhen  F'  is  DFR 


Sampling  with  replacement.  In  this  case 


^  nT 

G(T)  =  ^  = 


Ni(T) 


t  h 

".here  N.(T)  denotes  the  number  of  replacements  in  the  i  item 

n 

position  and  r  -  T  N.(T)  denotes  the  total  number  of  replace- 

1=1  '  r  1  1 

ments  in  f O,  T ]  .  Of  course  E - I  is  unbounded. 

E  N.(T) 


T 

However,  '"e  kno",-  that  E[  N.(T)]  for  all  T  ^  0  (Barb". 

A 

and  Proschan,  1964b).  Hence  this  again  indicates  that  0  (T)  v/ill 
tend  to  be  larger  in  the  IFR  case  than  in  the  exponential  case. 


4.  Bounds  on  time  to  r  failure.  Under  the  exponential  ass  ’mp- 

tion,  the  distribution  of  the  statistic  6  depends  only  on  r  and 

r,  n 

not  on  n  .  The  choice  of  n  in  this  case  is  usually  determined 
by  the  ratio 


E(X  ) 

_ Li.!L_ 

which  is  an  indirect  measure  of  the  expected  saving  in  time  due  to 
putting  more  than  r  items  on  test  but  terminating  at  the  r^^ 
failure  (Epstein,  1960a).  We  always  have 

<  1  . 

■ 


Since  the  bound  is  ritt.iineri  by  the  degen**  r.’ te  distributH)n  (  hicii  is 
the  limit  of  IFR  distributions),  this  is  not  a  useful  measur*'  if 
we  assume  only  IFR  .  Ho-’-ever,  "e  can  obtain  n on-trivial  bounris 


on  F(X  )  . 

r ,  n 


Assume  F  Is  IFR,  •'/ith  mean  1.  and  continuous.  We 

may  "  rite  Y  -  -locF^fX.)  where  Y.  is  the  i^^  order  statistic 
11  1 


-  \ 


,  and  is  ,i  con- 


in  a  sample  of  n  from  distribution  G(x)  -  1 

t  h 

vex  function  of  X.,  'here  X.  is  the  i  order  statistic  in  a 

i  1 

sample  of  n  from  F  .  By  Jensen's  inequality 

E(YJ  >  -log  r  fF(\X)1. 

so  that 

rr  F(x.)i  2  e'^'^'^i^  . 

If  b(x)  is  a  sharp  upper  bound  on  Fix)  ,  then  b  is  decreasing 
to  0  and 

bf  E(X^)1  >  e'^'^'^i^  . 

Hence  choosing  Xq  such  that 
b(xQ)  -  e  '  1 

vhere  of  course 


^  I  • 

we  have  j  =  l 

E(X_)  <  . 

Using  tabled  upper  bounds  on  F  given  one  or  tv  o  moments  of  F 
(Barlow  and  Marshall,  1963)  we  can  obtain  upper  bounds  on  E(X.)  . 
When  F  is  DFP  we  can,  in  a  similar  manner,  obtain  lover 


bounds  on  E(X.) 


using  lower  bounds  on  F 


.’J) 

If  'V  e  specify  the  first  nionient  of  F  ,  explicit  upper  Ivtunds 
can  be  given  on  F(X.)  '‘  hen  F  is  IFR  ,  as  sho  "  n  in 

T  heorem  4.1.  If  X  ^  X  ,  ^  .  .  .  -^  X  are  order  statistic  s  froni 

- 1  _  j  f-i 

an  IFR  random  variable  •.>ith  mean  0  and  Y ,  Y  ,  '  .  '  Y 

1  —  -  -  I, 

are  »he  order  statistics  from  G(x)  1  -  e  ^  .  then 


(a) 

r-F(Yj)  <  F{Xj)  - 

f- 

(b) 

F(X.)  - 

1  — 

f’  E(Y., 

1  ^  i  ^  n 

(ct 

0  <  F(X  )  < 
■"  n  — 

e  e(y  ) 

n 

(a)  and  (c)  are  sharp  and  (b)  is  non-trivial  though  not  sharp. 


Proof,  (a)  and  (c)  are  sho”  n  in  Barlo”  and  Proschan  (19b4a), 


Chapter  2.  Hence  e  need  only  oro/e  (b)  and  ve  may  assume  H  t  1 

First  let  us  verify  that  (b)  is  non-trivial.  Note  that  by  (a)  , 

EX  ,  <  EX  <  EY  ,  so  that  FY  is  a  trivial  upper  bound  for 
n-l—  n—  n  n 

EX  ,  .  Therefoie,  a  non-trivial  upper  bound  for  EX  ,  must  be 

n-l  n-l 


F  Y 


less  than  EY  ;  i.  e.  ,  .ve  must  sho  e  that 
n 


lUL -  -  FY 

-^^n-1  ” 


But  for  z  >  0  , 


1  -  e 


-  z 


1  -  e 

/.  *^1  ;  thus  letting  z  -  FY 


1-1 


1 


—  + - —  r-  +  .  .  .  "  e  conclude 

n  n  -  1  2 


FY 


n  -  1 

-FlYn-l 


1  -  e 

To  show  (b)  use  the  bound 
r(x)  ^  b(x)  = 


1 

—  f 

n 


1  i  1 

.rrr  ‘  *  7  ‘  ‘ 


RY 


X  ^  f’ 


VC  X 


X  >  f' 


.  X 


where  w  depends  on  x  and  satisfies  \  e  du  =  T  . 
(Barlow  and  Marshall,  19f>4).|| 

Sharp  bounds  will  be  derived  in  a  future  publication.  However 
these  are  not  as  convenient  as  the  bounds  of  Theorem  4.1. 

Bounds  on  expected  values  of  order  statistics  can  also  be 
given  in  terms  of  the  p^^'  percentile. 


Theorem  4.  Z.  Let  X,  <  X  ,  <  .  .  .  <  X  denote  the  order  statis- 
-  I  —  Z  —  —  n 

tics  from  F  ,  IFR  with  p*^  percentile  Then 


(4.1) 


E(X  )  <  maxy  F  ,  — E —  . 

J  -  I  P  -log  q 


1  1  \ 


(4.2) 


J-1 

e(x4  >  y  f  r 


j'-  ^ 
1=0 


(r  Jo 


X  log  q  X  log  q 


l-e  ‘e 


where  q  =  1  -  p  .  All  inequalities  are  sharp. 


Proof.  To  show  (4.1),  let 


0  <  X  <  A 


C^(x)  = 

9  e’‘p[  (  r  .  ^ 


Note  that  G^(A)  =  1  and  =  1  -  p  =  q  .  Since  log  F(x) 

is  concave,  there  exists  at  least  one  value  of  ^  >  0  such  that 


C^(x)  '»F(x)  forall  x^O  .  Thus  E(X .)  ^  sup  E(  Y .)  where 

■’  "o  <  A<c 

-  -  P 

Y.  is  the  J  order  statis  ic  from  G*  N aw 
J  ^ 

i-1 


E(Y, 


)  =  P  ^  P'ilG  (x)|‘lC  (x)r-‘dx 

“  *  I  rwnno’-j)  L  . 


'C^(x) 


To  find  the  maximizing 


A  ,  consider 


Since  C 

r(n 

TTITTVT 

where 

Hence 

-  ( 

For  j 

we  have 
case  at 

F.(Yj)  = 


1  - 


r(n  f  1) 

r(jjr(n  +  1  -  j) 


.  1 


t-^'V-t)"'-’  dt 


G^(A) 


*  rnyrurrrrj]  \  fr,(x)|J ‘(C(x)i 


n- 


X  -  ^ 


qexp[-p - ^  logq)  logq - 


:^(A)=0  . 


riE(Y  J 

-  -  A)  — w—— —  reduces  to 

P  r*  A 


t  1) 


OC 


g^(x)(x-5p)  dx  =  F(Y|)-«p 


g^  is  the  density  of  . 


?  -  A)  E{Y  )  =  A--E -  (if -  f  - i - p)  -  4 

^  p  -j  A  '  j  lo^g  q  '  n  n  -  j  f  1  p 

=  -(4  -  A)[  1  f  jJ—  {Li -  f - 1_)  1  . 

p  log  q  n  n - j  ^ 1 

such  that  1  4-  1 —  (  i  f  ....  f  - ! — ^  0,  we  have 

log  q  n  n-jfl 

<  0  .  For  1  such  that  1  *■  • — ^ —  (  i  f  .  .  .  .  +  — I — ^  )  >  0 
—  logqn  n-jfl  — 


St 


TrE(Y.)  >  0  .  Thus  E(Y.)  is  maximized  in  the  fir 

^  J  -  J 

A  =  0  and  at  A  =  4  i»^  the  second  case.  When  A  =  0, 

P 

—  -  (if - f  - )  .  When  A  =  4  ,  E(Y.)  =  4  • 

-  logq  n  n  -  j  f  1  p  J  p 


Let 


C(x)  =  / 


x  log 


for  0  <  x  <  4 


4  <  X  <  X 

P  - 


T o  show  (4.  Z). 


0 


for 


Then  G(x)  <  ^(x)  for  all  x  >  0  and  G  has  percentile 

^  .  Thus  E(Y  )  <  F(X  )  where  Y  is  the  order  statistic 

P  J  ~  J  J 

from  G-.  But  , 

r'-V  ' 

E(Y  )  =  I  )  "  (G(x)l'  dx  . 

'  ito 

Using  the  above  definition  of  G  ,  we  obtain  (4.  Z). 


5.  Further  results  on  order  statistics  and  spacmgs.  In  this  section 
we  pi  esent  some  results  of  theoretical  interest  concerning  order 
statistics  and  their  spacings  from  and  IFR  (DFR)  distri¬ 

butions.  Many  of  the  results  hold  without  the  restriction  F(0  )  =  0  , 
First  we  present  some  total  positivity  properties  of  the  order 
statistics.  A  function  K(x,  y)  of  two  real  variables  is  said  to  be 
totally  positive  of  order  r  (TP^)  if  for  all  1  <  m  <  r; 

X,  <  X  ,  <  .  .  .  .  X  and  y,  <  y y  we  have  the 
1—  I  —  —  m  ’  \  ’  I  -  —  'm 

determinant  inequalities 


K(x^,  y^) 


m 

i.  J=1 


>  0  . 


The  following  lemma  is  of  use  later  on,  as  well  as  of  interest  in  its 
own  right. 


Lemma  5.1.  Let  F  be  a  distribution  having  density  f  with 
f(x)  not  necessarily  0  for  x  <  0  .  Let  f^(x)  be  the  density 

of  the  i^^  order  statistic  in  a  sample  of  size  n  .  Then  f^(x)  is 

TP  in  i,  X  where  i  =  1,  . n  and  -  oo  <  x  <  «  . 

D 

Proof. 


(5.1) 


Since  -  — ' —  is  T  F’  in  1  and  x  ,  whi-n  i  1 ,  Z ,  .  .  .  ,  n 

r(x)^ 

1 1 

and  -  00  ^  x  *'  at  ,  the  conclusion  follows.  i  I 

We  may  obtain  a  similar  result  concerning  the  right  hand  tail 
of  the  distribution  of  an  order  statistic. 

Lemma  ^.2..  Let  F  be  <  ny  distribution  with  F(x)  not  neces¬ 
sarily  0  for  X  <  0  ,  F.  the  corresponding  distribution  of  the 

i^^  order  statistic.  Then  F  (x)  is  TP  in  i,  x  ,  where 

1  ur 


i  =  1,  . n  and  -  -r  ■c  x  <  or  . 

Proof.  .  , 

-  i- 1 

(S.2)  F(x)  =  ^  F^lx)  f"*J(x)  . 

j=0  ^ 


r  F(x) 

fsjow  J  — ■ — -  ,  is  TP  in  x,  j  .  Therefore 
^  F(x)  '  * 

00 

F.  (x)  =  y  r"(x)(-^^^’'  H(i  -  1  .  j)  IS  also 

‘  j_Q^  ^  ^  r(x)  ^ 

TP  in  i  and  x  .  where  H{k)  -  1  for  k  >  0  ,  0  otherwise.  I  ' 

or  — 

In  Barlow,  Marshall,  Proschan  (1963),  it  is  shown  that  the 
order  statistics  from  an  IFR  distribution  themselves  have  an  IFR 
distribution.  The  next  lemma  shows  a  similar  preservation  of  the 
PF^  property. 

Lemma  5.3.  Suppose  the  underlying  density  f  is  FF^  ,  with 
f(x)  not  necessarily  3  for  x  <  0  .  Then  the  density  f.  of  the 
i^^  order  statistic  is  also  PF^  for  fixed  i  =  1,  Z, . n 
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Proof.  It  is  easy  to  verify  that  when  f  is  ^^2  »  ^  and 

F  .  Thus  log  f,  log  F,  and  log  F  are  concave.  It  follows  from 
(5.1)  that  log  f^  is  concave,  or  equivalently,  f^  is  for 

fixed  i  =  1,  2 . n  .  i  j 

Next  we  obtt  .i  some  corrparisons  between  the  order  statistics 
of  an  IFR  (DFR)  distribution  and  the  corresponding  order  statis¬ 
tics  of  an  exponential  distribution. 

Theorem  5.4.  Let  X,  <  X,  <....<X  ,n>2,be  order  sta- 

-  1  —  2  —  —  n  — 

tistics  from  F  ,  an  IFR  (DFR)  distribution  with  mean  6  ,  but 
t 

^  1  -  e  .  Let  Y,  <  Y  ,  <  .  .  .  .  <  Y  be  order  statistics  from 

1  —  2  —  —  n 

t 

G(t)  =  1  -  e  .  Then 

(a)  EXj  -  EY^  has  at  most  one  change  of  sign  as  j  goes 
from  1  to  n  .  Moreover  if  one  change  of  sign  does  occur,  then 
EXj  -  EYj  goes  from  positive  (negative)  to  negative  (positive)  values. 

(b)  If  6  =  6',  then  one  change  of  sign  does  occur. 

Proof. 

(a)  Assume  F  is  a  continuous  IFR  distribution  with 

mean  6  and  G  is  exponential  with  mean  6*  .  We  have  seen  in 

Y 

Section  1  that  if  Y  has  distribution  G  ,  then  X  =  <>(-^)  has  dis¬ 
tribution  F  ,  where  4)”^(x)  =  -logF(x)  ,  a  convex  increasing 
function  which  is  0  for  x  =  0  . 

Thus 

EXj  lo  *  y}gi(y)dy  • 


where  j’.  is  the  density  of  the  i  order  statistic  from  tht  ex¬ 
ponential  distribution  G  .  By  Lemma  S.l,  R-ly)  ^  ^  cr,  '''  ' 

and  V  .  Also  -  y  changes  sign  at  most  once,  and  if  once, 

from  positive  to  negative  values.  Py  the  variation  diminishing 
property  of  totally  positive  functions  (Karlin,  19fi4,  p.  LI),  F'X.  -  P. 
also  changes  sign  at  most  once  ,  and  from  positive  to  neg<iti\'e  va  lu  es, 
if  at  all. 

n  n 

(b)  If  (<  =  O' ,  then  1  EX.  =  fF  =  h'  r  T.  PlY.  .  Hence  P'.X,  -  PV 

.  ,  1  1  1  1  ' 

1=1  1=1 

must  change  sign  at  least  once  or  be  identir  ally  0  fr.. 

i  =  1,  2 . .  n  .  Now  since  F  G  ,  F  cannot  agree  with  G 

on  an  interval.  Hence  by  Corollary  ‘1.  10,  Chapter  2  of  Barlow  and 

Proschan  (1964),  EX,  >  EY,  and  EX  <EY  .  Thus  EX  -  EY 

11  n  n  '.1 

is  not  identically  0  for  i  =  1,  2,  .  .  .  .  ,  n.  Hence  EX.  • 
chai.ges  sign  exactly  once. 

If  F  is  IFR  but  not  continuous,  we  niay  obtain  the  same 

result  b^  using  continuous  IFR  approximations. 

Finally,  a  similar  argument  holds  if  F  is  DFR.  Ij 

tually,  under  the  same  hypothesis  we  may  prove  a  stronger 

version  of  (a)  in  which  P'.X.  -  EY.  is  replaced  by  FiiX^  -  P'Y^  , 

o  '>  0  .  Jf  instead  of  assuming  EX  =  EY  in  (b),  we  assume 

EX^  “  EY^  ,  then  we  may  show  that  one  change  of  sign  of  EX^  -  EY^ 

does  occur.  We  omit  tne  details. 

We  may  obtain  further  consequences  of  Theorem  S.4  using  the 

notion  of  major ization .  A  vector  a  =  (aj,a^ . a^)  m.ajorizes 

a  vector  b  =  (b,,b  . b  )  (written  a  >  b)  if 

—  i  c  n  —  — 


k  k 

a,  >  a  >  .  .  .  .  >  a  ,  b  >  b  >  .  .  .  .  >  b  ,  1  a  >  ^  b. , 

I  —  c  —  —  n  1—  Z—  —  n 

i  -1  1=1 

n  n 

k  =  1,  Z,  .  .  .  .  ,  n- 1 ,  and  Z  a  =  2.  b.  .  See  Hardy,  Littlewt)od, 

.  ,  i  .  ,  1 

1=1  1=1 

Polya  (19‘^t^),  Chapter  II. 


Theorem  S.  5.  Let  X,  X 

- 1  —  —  n 


be  order  statistics  from  F 


an  IFR  (DFR)  distribution  with  mean  (>  ,  Y  Y  be 

1  —  —  n 

t 

order  statistics  from  G(t)  -  1  -  e^  .  Then  (EY  ,  EY  . . EY.) 

n  n- 1  1 


>(-)(EX  EX 

n  n  - 1 


.EXj)  . 


Proof.  Let  F'  be  IFR  .  From  Theorem  5.4  we  know 

EY^  -  EX^  has  one  change  of  sign,  from  plus  to  minus  as 

n  n 

i  goes  from  1  to  n  .  We  also  know  E  EY  .  ,  =  n&  =  1  EX  .  , 
“  .  ,  n-i4l  .  ,  n-i>l 

1=1  1=1 

J  J 

Thus  r  EY  .  ,  >  I  EX  .  j  =  1.2 . n  .  Finally,  EY  .  , 

.  ,  n-ifl— .  ,  n-ifl  '  n-i4l 

1=1  1=1 


and  EX  .  ,  are  decreasing  in  i  .  Thus  the  conclusion  follows, 
n-i+1  ^ 

A  similar  argument  holds  if  F  is  DFR  . 


Using  Karamata's  Theorem  we  obtain  Theorem  5.6  below. 
Karamata's  Theorem  states  that  if  0  is  continuous  and  convex  and 
a  ^  b  ,  then 

—  —  ’  n  n 

(5.  3)  >  )  t'(b.)  . 

i  =J  I'Tl 

See  Haroy,  Llttlewood,  Pt^lya  (195Z),  p.  89. 


Theorem  5.  6.  Let  il  be  continuous  and  convex, 


Y,,  .  .  .  .  ,  Y  as  in  Theorem  5.  5.  Then 
1  n 

n  n 


^0(EY.) 

i-1 


>20(EX.)  . 

i=l 


X 


( 

n 


Proof.  1  et 

P 

be 

IPH  . 

■'H 

I  eor  ( 

•  I  n  ‘1 

.  S, 

.  FYj) 

>  (FX  ,  .  .  . 

n 

.  ,  FX 

i>  ■ 

!  i  enc  e 

by 

Ka  r  <1  rna  ta '  s 

The 

■  or  ern , 

the 

c  one lu - 

sion  follows 

1.  A 

s  1  m  1 

lar  argurne 

nt  holo: 

s  if 

P' 

is  DP 

■R 

.  1 

Using 

rhe> 

)r  ern 

h  .  S  xt 

•  obt 

a  1 

Theorem  S. 

7.  I 

.et 

c  ■  c 

X,, 

. X 

Y . . Y 

'  ' 

1  - 

L-  — 

— 

r'i 

1 

n 

1 

n 

n 

as  in  Theor 

e  m  S , 

.  S. 

1  hen 

1  c 

F  Y 

c 

,  FX 

1  1  - 

1 

i 

i 

-1 

1  1 

Proof.  Let 

P' 

be 

IP’R  . 

Defining 

d  = 

1 

^  |- 
n  -  1  •  1 

FX 

r'.  -  1  » 1  ’ 

write 

n 


n 


^  c  .  ,(F.Y  ,  -  F.X  .  ^  (c  -  I  ,)ci,  .  (i  ,  -  c  i(d,  ^  d,) 

/_  n-14-r  n-i»l  n-it-r  '  n  n-l  1  '  n-1  n-Z  '  1  I 

1=1 

■  <S,.Z  ■  S,-3>'''l  •  'I,  ■  d,»  -  •  (c,  •  ..  .  -d,,.,) 


t  .....  dj 


Since  c,  -  c.  ,  0  ,  i  =  1 . n-l,  d  •  .  .  .  f  d  >0  , 

1  1  *■  1  —  1  I  - 

j  =  1 . .  n-l,  and  d^  *  .  .  .  *  d^  -  0,  we  rone  lude  that 

n 

c  .  ,  (I-.Y  .  ,  -  FX  .)  >  0  . 

/  n-iM  n-i*l  n-i‘i  — 

T-i 

A  similar  argument  holds  v«.’hen  F  is  I^FR.  I| 

Finally  we  sumtiari/.e  some  results  concerning  the  covariance 
of  order  statistics  cjbtained  by  Tukey  (IdSb).  He  -nows  ttiat  if  F 
is  IFR  ,  thenfor  h^'j-  k  .  cov  (X,.  XJ  '■  cov(X  .  .X  )  ,  where 

—  ■  k  h  J  b 

X  is  the  order  statistic  from  P'  .  He  furthc,'  shows  that  if 

,1 

F  satisfies  both 

(a)  log  F  is  concave  (i.  e.  ,  P'  is  IP"R),  and 

(b)  log  F  IS  concave, 
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then 

(1)  the  covariance  of  any  two  order  statistics  is  less  than  the 
variance  of  either,  and 

(2)  the  covariance  between  order  statistics  X.,X,  is  mono- 

J 

tone  in  j  and  k  separately,  decreasing  as  j  and  k  separate 
from  one  another. 

Note  that  if  f  is  ^  satisfies  both  (a)  and  (b)  above. 

Next  we  derive  properties  of  the  spacings  X^,  X^  -  X^ . 

X  -  X  ,  from  PF,  and  IFR  (DFR)  distributions  similar  to 
n  n-1  2 

those  of  the  order  statistics  X,  <  X,  <  .  .  .  <  .X  obtained  above. 

i  —  Z  —  —  n 

We  first  consider  total  positivity  properties. 

Theorem  5.8.  Let  f  be  P^2  neceasarily  0 

for  X  <  0  .  Then  h.  ,  the  density  of  X.  -  X.  ,  ,  is  PF  for 

fixed  i  =  2,  3,  .  .  .  ,  n.  If  we  assume  further  that  f(x)  =  0  for 
X  <  0  ,,  then  h^  is  Pf^2  ’  density  ol  X^  . 

Proof. 

(5.4)  ^i^^^  ~  (l-2)l(n-l)  *  ^(u)f(u)f(u  f  x)  F"  ^u  +  x)  du 

for  i  =  2,3 . n  . 

Since  f  is  PF^  .  so  is  r(u)  =  F^  ^(-u)f(-u)  , 

s(u)  =  f(u)  F^  ^^\u)  .  Hence  so  is 

N*-’  =  (i-DiVn-ia)!  ^  >‘i  0“ 

for  fixed  i  =  2,3 . n  . 

Assuming  f(x)  =  0  for  x  <  0  ,  wm  see  that  h^  is  ^F^,  from 
hj{x)  =  n  f(x)  F.  "  \x)  .  II 


(5.5) 
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Theorem  5.9.  If  P"  is  DP'R  ,  then  H  is  DP'l^  for  fixed 
-  1 

i  =  1 ,  ,  .  .  .  ,  n 

Proof.  Since  I)F!f  is  preserved  under  convex  combinations 
(Barlow,  Marshall,  P^roscluin,  19b  3,  p.  381)  we  see  from  the  repre¬ 
sentation 

(5.6)  \  f(u)  (x)  du. 

where  .  I^P  if  m  x  for  fixed  u  ,  that  !1.  is 

DFR  for  fixed  i  =1,^3 . ri*  l| 

Theorem  5.10.  Let  F  he  DFR  with  P'(x)'  1  for  all  x  >  0  . 
Then  H.(x)  is  where  i  =  <3,  3,  .  .  .  ,  n  and  x  >  0 

Proof.  F^  ^(u)  f(u)  is  TP  ,  in  i  ^  ,  3,  .  .  .  ,  n  and  u  >  0. 

^^\ut-x)  is  TP,  in  i,u,in  i,x,andinu,x(u>0,  x>0). 
Thus  by  a  theorem  in  the  book  by  Karlin  (in  process) 

(5.  6)  rr.(x)  =  ({  I)'i‘(n  IV  ^  ^ 

is  in  i,  X  ,  where  i  =  Z,  3 . n  and  x  >  0  . 

❖  <  _ 

Theorem  (Karlin)  let  \,  x,  ^  traverse  linear  sets  ,  X  ,  ^ 

respectively.  Suppose  h(\,  x)  -  ^  f(\,  x,  <  )  g(\,  ^')  d|jL(^)  is  well  de- 

fined  on  A  x  X  ,  where  p  is  a  tr  -finite  measure,  and 

(i)  i(\,x,t)  >  0  for  all  \  in  \  ,  x  in  X  ,  and  ^  in 

^  ,  and  f  is  TP^  for  each  pair  of  variables  when  the  third 
variable  is  held  fixed. 

(ii)  g(X.?)  IS  TP^  . 

Z  ' 


Then  h  is  TP 


il 


Next  Wf  present  sorin  in<i  t<<r  i /.<it  ion  properties  of  the  nor¬ 
malized  sp.K  inps  (n  -  1  •  1)(.\  -  \  j)  ,1  1 ,  . n  ,  similar  to 

ih.'  .Se  de\'elo[je(l  above  for  the  order  st.itistn  s. 

Theorem  S.ll.  Let  ^  X  ,  .  .  X  be  th--  order  statistics 

from  F  ,  an  Ih'K  distribution  .vith  mean  \ ,  Y  '  .  .  ■  V 

^  1  ~  —  n 

'  r< 

the  ord<  r  statistic  s  from  (j(0  1  -  e  .  Then 

FnX,,  !.(n  -  1)(X,  -  X  ) . !:(X  -  X  ) 

1  ..1  n  n  -  1 

>(-)  FnY  F(n  -  1)(Y  ,  -  Y.) . F(Y  -  Y  , )  . 

I  ^  i  n  n  - 1 


Proof.  Let  h'  be  IFi^  .  Bv  Theorem  1,  i'(n  -  i  *  1)(X^  -  X^  j) 

IS  decreasing  in  i  ,  It  is  alsc:)  easy  to  verify  th.it 

-  F(n  -  1)(Y  ,  -  Y.)  =  .  .  .  =  f;(Y  -  Y  .  Since 

z  1  n  n  - 1 

n  n 

F(n  -  1  ‘  1)(X^  -  X^_j^  --  b  F  (n  -  i  ^  1)(Y^  -  Y^_j)  , 

1 1  T- 1 


it  follows  that 


2^(0 

1=1 


-  1  *  l)(X^  -  >  F.(n  -  1  .  1)(Y,  -  Y^_j) 


for  j  =  I,  Z  ,  ,  n-1  .  Thus  the  c'onrlusion  follocvs. 

A  similar  argument  holds  if  F  is  DP'R.  | 

For  normalized  spacings,  the  analogue  of  Theorem  S.7  is 

Theorem  S.IZ.  Let  c,''...  >c  ,  X,.X,...,X  ,  Y,,Y, . Y 

-  1-  —  n!z  nl.;  n 

as  in  Theorem  S.  11.  Then 

n  n 

y  CjE(n-i  M)(X^  -  X^_j)  2  ")  CjE(n  -  1  sl)(Y^-  Y^_y  . 

1  =  1  (•')i=l 


Proof.  The  proof  parallels  that  of  I'hcorem  5.7. 


We  immediately  obtain: 

CorollaryS.H.  Let  X, . X.Y, . Y  be  as  in  Theorem  5. 11. 

— - - - -  1  n  1  n 

Then  for  I  r  <'  n 

r  r 

^  E(n-i  M){X.  -  X^_^)  >  ^  E(n  -  i  f  1)(Y.  -  Y^_j)  . 
i=l  i“l 

Proof.  Choose  t,  =  c,=...  -  c  =1,  c  ,=c  ,-...  =  c  -0  , 

-  1  Z  r  ’  r  >1  r  n 

so  that  c,  >  c  >  .  .  .  >  c  ..  The  result  follows  from  Theorem  5.  11. 

1  —  Z  -  —  n 
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on  tha  roport,  aaa  data  of  publication 


lapoood  by  aacurity  claaalflcation,  uolng  atandard  statononto 
auch  aa: 

(1)  “Qaallflod  raqaaataro  aay  obtain  copiaa  of  thla 
roport  hoai  DDC*' 

(2)  '*Porol0i  anaouncaaMat  aad  dIoaoaUnatlon  of  thla 
roport  by  DDC  la  aot  anthortaod  " 


(3)  “U.  E  QovarBMaal  aeoaclao  aay  obtain  copiaa  of 

thla  report  diractly  froai  DDC.  Other  quail  od  DDC 
uaara  ahall  raquaot  through 


(4)  **U.  E  adlltary  agonclao  any  obtain  copiaa  of  this 

roport  diractly  front  DDC  Othor  qualifiod  uaara 
ahall  raquaat  through 


(5)  "AU  dlatrlbution  of  thla  report  ia  eontrolloE  Qual- 
ifiad  DDC  aaaro  ahaU  raquaat  through 


If  tha  lapaet  haa  boaa  faralohad  to  tha  OfSco  of  Technical 
Sarvlcaa,  PigartaiaBt  of  Commmem,  for  aala  to  tha  public,  ladl- 
eata  thla  fact  aad  aatar  tha  prica.  If  known 

IL  SUPPLEMENTARY  NOTBE  Uaa  for  addlUoaal  oaplano. 


IX  SPONSORING  MILITARY  ACTIVITY:  Eatar  the  name  of 
tiM  dapartnoalal  projact  ofSca  or  laboratory  apoaaorlog  (pay 
fas  lot)  tha  raaaarch  and  davalepntaat  lacludo  ad*aan 


7a.  TOTAL  NUMBER  OP  PAGER  TIm  total  p^  count 
ahould  follow  aaeaial  pagtaMloa  procadurao,  Ln,  aatar  tha 
number  of  pagaa  rnetahilag  lafotiaatlon 

7b.  NUIBBR  OP  REPBRENCBR  Eatar  tha  total  auahor  a 
rofarancaa  cltod  la  tha  report. 

So.  CONTRACT  OR  GRANT  NUIBBR:  If  appr  opr  lata,  aatar 
tha  applicabla  Maher  of  tha  conhact  or  graat  aadar  trtticb 
tha  roport  araa  ealttan 

Sb.  R,  R  Sd.  PROJECT  NUMBER:  Eatar  tha  agprapriata 
military  iopaNmaat  Idaatlflr  atloa,  auch  aa  prajM  aambar, 
aufeproiact  aaaRar,  ayatan  aaRbara,  tank  naaRar,  ate. 

9a.  ORIGINATOR'S  REPORT  HUIBERfS):  BMar  the  ofE* 
clal  report  aaabar  by  oddeh  tha  decanaat  erUl  ba  Idaatlflad 
aad  contraUad  by  tha  orlglaatlag  activity.  Thla  naRbar  anu 
ba  unlqua  ta  tUa  lapart. 

9».  OTHER  REPORT  NUtBERfS):  U  tha  roport  haa  baaa 
aoalgnad  aay  othor  report  aMhara  fafthar  by  tha  arfglaator 
or  by  (ha  epaaaerX  alaa  aatar  thla  nuMbaKa). 

10.  AVAlLABaJTT/LBRTATKNI  NOTICBR  Eatar  any  U 
itatloaa  oa  fthar  diaaamlaatina  af  tha  rapart,  oSmt  thaa  tl 


19.  ABSTRACT:  Baler  aa  abatiact  glviag  a  brief  and  factual 
oaiMary  of  tha  docaaaat  ladlcatlvo  of  the  report  avaa  though 
It  may  alao  appear  alaawhaia  la  the  body  of  tha  tachaical  ra¬ 
part.  If  additloaal  apace  ia  raquirad.  a  coatiauation  ahaat  ahall ' 
ba  attached. 

It  la  hlBIy  daairablo  that  the  abatract  of  claaalfiad  roporta 
ba  aaelaaelAad.  Bach  parn^ph  of  Sm  obatract  ahall  aad  with 
aa  Indication  of  the  miUtary  aacurity  claaalflcation  of  tha  la- 
fonaatlaa  ia  the  paiagyaph,  lapraaoatad  aa  (Tt).  rs>.  fc>.  ar  (V). 

There  la  no  lialtatlon  on  tha  laagth  of  the  abatract.  How- 
aver,  Sm  auggaatad  laagth  la  frooi  ISO  to  223  wordo. 


14.  BBT  BOROS:  Bay  ararda  are  rachalcally  ataaalagftil  taima 
ar  ahort  phiaaaa  that  charactariaa  a  (apart  aad  may  ba  uaad  aa 
ladaa  aatrtaa  far  eatalaglag  Sm  rapart  Bay  worda  muat  ba 
aalactad  aa  SMt  aa  aacurity  claaalficatlea  ia  raqvUfod.  Idantl- 
flara.  aaeh  aa  aqulpMat  ■ado I  daalgiatloa.  trade  nama,  eUUury 
ptaiact  eada  aaaa,  arovnphlc  location,  UMy  bo  good  aa  hay 
trarda  bat  trlU  ba  faUaerad  by  an  ladicatloa  of  tachaical  coa- 
taat.  Tha  aaolgamiat  af  tlaka.  lalaa,  aad  tral^a  la  opUonal. 
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